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ABSTRACT 

Using the Subaru 8.2m Telescope with an IRCS Echelle spectrograph, we obtained high-resolution 
(R =10,000) near-infrared (1.01-1.38/xm) spectra of images A and B of the gravitationally lensed QSO 
B1422+231 (z = 3.628) consisting of four known lensed images. We detected Mg II absorption lines 
at z = 3.54, which show a large variance of column densities (^0.3 dex) and velocities (~10 km s -1 ) 
between the sightlines A and B with a projected separation of only 8Ah^Q pc at the redshift. This is 
the smallest spatial structure of the high-z gas clouds ever detected after Rauch et al. found a 20-pc 
scale structure for the same z = 3.54 absorption system using optical spectra of images A and C. 
The observed systematic variances imply that the system is an expanding shell as originally suggested 
by Rauch et al. By combining the data for three sightlines, we managed to constrain the radius and 
expansion velocity of the shell (^50-100 pc, 130 km s _1 ), concluding that the shell is truly a supernova 
remnant (SNR) rather than other types of shell objects, such as a giant H II region. We also detected 
strong Fe II absorption lines for this system, but with much broader Doppler width than that of 
a-element lines. We suggest that this Fe II absorption line originates in a localized Fe II-rich gas cloud 
that is not completely mixed with plowed ambient interstellar gas clouds showing other a-element 
low-ion absorption lines. Along with the Fe richness, we conclude that the SNR is produced by an 
SNIa explosion. 

Subject headings: galaxies: abundances-gravitational lensing: strong-intergalactic medium-ISM: su- 
pernova remnants-quasars: absorption lines-quasars: individual (B1422+231) 



1. INTRODUCTION 

Gravitationally lensed QSOs provide us precious op- 
portunities to study spatial structures of high-z gas 
clouds which int ersect the multiple sight li nes toward 
the QSO (e.g.. iWevmann fe Foltzl [19831 iFoltz et alj 



1984; 
2002: 



Smette et aL. .. . 
_ _ .Kobavashi et all I2002T . 
Ellison et al.l l2004t lLopez et al 



19951: IRauch et all I1999L~ 120013111: 
Churchill et al.lT2 003a: 
_ ; _ 120051 : iMonier et al.l 

20091 ). By comparing the profiles of absorption lines be- 
tween multiple sightlines, we can study the density and 
velocity gradients of intersecting gas clouds directly even 
at high redshift. Such spatial properties of absorbing gas 
at high redshift may provide a key to understand the 
galaxy formation processes. 

The statistical studies of absorption line systems 
with the gravitationally lensed QSOs have revealed 
that the Mg II systems, which trace low-ionization 
systems are relatively small (< a few hundred pc) 



and have clumpy s patial structures (jRauch et al.l 120021 
lEllison et al.1 H004 1 ) , in contrast to the C IV systems 
which are typically a few kpc in size and have fewer struc- 
tures ([Rauch et al.ll2001atfEli"ison et al.ll2004t l. Studies of 
single line-of-sight observations and C LOUDY photoion- 
ization modeling ([Ferland et al.lli~998[ ) also suggest that 
the Fe -rich low-ionization systems should be as small as 
10 pc (|Rigbv et al.l [20021 liCravanan et al.|[2008h . Thus, 
low-ionization systems appear to have small clumpy spa- 
tial structures that directly relate to star forming activ- 
ities, such as giant molecular clouds. It is important to 
investigate their spatial properties with multiple sight- 
lines of gravitationally lensed QSOs, since only the grav- 
itational lens can directly reveal fundamental parameters 
such as the size and kinematics of gas clouds. 

The quadruple gra vitationally lensed QSO B1422+231 
(jPatnaik et al.lll992l ) is one of the best objects for such a 
study because of the relatively small separations among 
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Fig. 1.— HST F555W image of B1422+231 from CASTLES 
homepage (http://www.cfa.harvard.edu/castles). North is up, east 
is to the lett. tour white point sources show gravitationally lensed 
images A, B, C and D. The extended object (G) to the northwest 
of image D shows the lensing galaxy of this system. We obtained 
spectra of the two brightest and closest components, A and B. 

multiple sightlines due to the large distance between the 
lense galaxy (z = 0.339; IKundic et all 119971 iTonrvl 
IMl) and the QSO (z = 3.628). This QSO has been 
observed frequently as one of th e most luminous grav- 
itationally lensed high-z QSOs (iBechtold fc Yed 119951: 
Songaila fe Cowiel 119961 iPetrv et all 119981: IRa uch et al l 
19991 l2001allbD . IRauch et al.l (|1999l , hereafter RSB99) 
obtained high- resolution optical spectr a of images A and 
C of B1422+231 using Keck HIRES (jVogt et al.lll99l 
and detected C IV, Si IV, C II, Si II, and O I absorp- 
tion lines fro m the z = 3.54 absorption system originally 
identified by iSongaila fe Cowiel (|1996l ) with strong Lya 
and Ly/3 absorption lines. The projected separation be- 
tween A and C sightlines is only 22.2/i^ 1 pc at the red- 
shift. RSB99 suggested that the absorption system is an 
expanding shell of mass ejection or supernova remnants 
(SNRs) by analyzing the differences of the absorption 
profiles of lower ionization species between the images A 
and C. 

In this paper, we report the results of the Subaru near- 
infrared spectroscopic observations of images A and B of 
B1422+231. Previous studies of low ionization gas with 
Mg II QSO absorption lines were limited to the optical 
wavelength range and thus to redshifts < 2.5, beyond 
which this transition moves into the infrared waveband. 
Thanks to the advent of highly sensitive near-infrared 
high-resolution spectroscopy with 8-10 m class telescope, 
it is possible now to study this unexplored redshift range. 
We are conducting a systematic high-resolution spectro- 
scopic survey of high-z Mg II systems with the Subaru 
IRCS Echelle spectrograph. We observed B1422+231 as 
one of the initial targets, and detected Mg II doublet and 
Fe II absorption lines for the z = 3.54 system. Moreover, 
we succeeded in spatially resolving spectra of images A 
and B owing to the high spatial resolution in the infrared 
and the very good seeing of the Subaru Telescope site. 
This paper is composed as following. The details of our 
observations are summarized in §2. The data reduction 
and calibration of spectra are described in §3. In §4, we 
show the properties of detected Mg II and Fe II absorp- 
tion lines. We interpret the properties as signatures of a 



Type la supernova (SN la) remnant, which is discussed 
in §5 in detail. §6 is the summary of this paper. We 
adopt a standard cosmology with = 0.7, VL m = 0.3, 
Ofe = and H = 70h^o — 70 km s _1 Mpc -1 throughout 
this paper. 

2. OBSERVATION 

We observed images A and B of B 1422+231 (Fig- 
ure Q} using Subaru 8.2m Teles cope (|Ive et al.l 120041) 
with IRCS3 Eche l le spe ctrograph (jTokunaga et al.lll998t 
iKobavashi eFaLI 12000) in Y (1.01-1. 19/zm) and J 
(1.16— 1.38^m) bands. The data were obtained on 2003 
February 13 and 2002 April 28 for Y and J bands, respec- 
tively. The weather condition for both observing runs 

was photometric and the seeing was good (~ .5) dur- 
ing observations for both bands. The integration time 
per frame was 600 sec, and the total integration time 
was 9000 s and 9600 s for Y and J bands, respectively. 
The slit widths were .60 and .30, and corresponding 
spectral resolutions (R — A/AA) were 5,000 and 10,000 
for Y and J bands, respectively, while the slit length 
was 3 .47 for both observations. Although we used the 
Subaru Adaptive O ptics 36-elements (AO 36) systerrQ 
(jTakami et al.| [2004) for the V-band observation, the im- 
provement of image quality was not good enough that we 
used the wider slit. The slit position angle (P-A.) was set 
at P.A. = 39°2 to put both A and B images on the slit 

simultaneously. The telescope was nodded by 1 .5 arcsec 
along the A— B direction as shown in Figure [H between 
alternating frames to offset the sky emission in the sub- 
sequent data reduction. We call the frames for positions 
"a" and "b" as shown in Figure [2] We also observed a 
standard star, 10 Boo, for flux calibration and removal 
of telluric absorption lines for both bands. The airmass 
of the target was distributed in a wide range from 1.0 to 
1.5 for both bands, while that of the standard star was 
about 1.0 and 1.5 for Y and J bands, respectively. 

3. DATA ANALYSIS 

3.1. Reduction 

We used IRAE0 for data reduction following standard 
procedures. First, we subtracted "b" frames from cor- 
responding "a" frames to cancel out the dark counts 
and the sky OH emission lines. All the subtracted 
frames were combined after the flat-fielding correction 
and hot/bad pixel correction. Then, we extracted two- 
dimensional (2D) spectra with the spatial axis along the 
slit and the dispersion axis perpendicular to the slit for 
each cross-disperser order using IRAF task "apall" in the 
"echelle" package. The 2D spectra of each order were 
combined after applying the wavelength calibration with 
Argon lamp spectra, which were obtained after the ob- 
servation. 

3.2. Deconvolution of A and B Spectra 

1 IRCS was mounted at the Cassegrain focus at that time. Now 
it is located at the Nasmyth focus. 

2 The A036 system is now upgraded to A0188 system with 
much improved image correction capability. 

3 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Fig. 2. — Upper panels: the slitviewer X-band images obtained during the echelle observation. These are raw images with an integration 

time of 120 s without corrections for hot/bad pixels. The pixel scale is .075 pixel -1 . The horizontal black thick bar near the center of 
the images is the shadow of the slit, and four point sources which are seen in white are the images A, B, C and D of B1422+231. The left 
and right panels correspond to the dither patterns "a" and "b" (see the text for detail) . The contamination of image C in the slit is slightly 
seen in this image. Lower panels: the spatial profiles of the spectra along the slit. Points show the observed data. Dotted lines show fitted 
profiles for each image A, B, and C, while thick lines show the combined profile of the three images. To show the spatial profile clearly, 
the data are averaged in the dispersion axis over one cross-disperser order (order 44 of J band). The hollows at both edges show the areas 
beyond the end of the slit, which are not used for the profile fitting. 



The lower panel of Figure [5] shows the spatial profile of 
the 2D spectra of B1422+231 along the slit. The images 
A and B are almost resolved, but the overlap is not neg- 
ligible. Image C, which is not a target of this study, is 
contaminated slightly in the spatial profiles. To obtain 
precise one-dimensional (ID) spectra of images A and B 
from the 2D spectra, we fitted the observed spatial profile 
with the point spread functions (PSFs) of the three im- 
ages pixel by pixe l in the dispersion a xis. This method is 
similar to that of lLopez et al.l (|2005l ). who also observed 
the multiple images of a gravitational-lensed QSO in a 
slit simultaneously. 

Because the tail of PSF was found to decline slower 
than that of a Gaussian function, we assumed the fol- 
lowing double Gaussians as the form of a PSF: 



fi(x,X) 



ai(A) 



exp 



exp 



{x-hY 



(i) 

where i —A, B, C are the indices for the lensed images, x 
is the coordinate of pixels in the spatial axis, dj(A) is the 
peak value of each image at wavelength A, hi is the center 
position of each image, and a s and ad are the width of 
each Gaussian function (a s < ad ; the subscripts "s" and 
"d" mean sharp and diffuse, respectively) . This function 
was found to fit the observed profile precisely as shown 
in the lower panel of Figure [2] 

The parameters bi, a s , and ad have to be fixed for 
linear fitting, which is necessary to converge the fit for 
observed profiles with low signal to noise ratio (S/N). 
These parameters are successfully determined from the 
fit of the high S/N spatial profile made by summing the 
spatial profiles for 50 pixels along the dispersion axis. 



The "fit" command of gnuplot was used for the fitting 
procedure. Though this program is not designed for nu- 
merical fitting, it is good enough to conduct the linear 
fitting. The flux density count at wavelength A for im- 
age i was determined by the integral of the PSF fi(x, A) 
along the spatial axis, which can be analytically calcu- 
lated as ^/Trai(X)(a s +ad)/2 for double Gaussians. When 
all the fitting parameters are determined for the entire 
wavelength range, the final ID spectrum of each image is 
constructed by plotting the flux density count along the 
dispersion axis. 

3.3. ID Spectra 

Many telluric absorptions of e.g., water (H2O) and 
oxygen (O2) plus emission lines of OH appear on the 
near-infrared spectra for ground-based observations. Al- 
though the emission lines are removed by the substrac- 
tion of images of alternative pointing ( §3.ip . the absorp- 
tion lines are still superposed on the extracted spectra 
of the object. For both bands, we removed the telluric 
absorption lines by dividing the object spectra by that of 
a AO standard star, 10 Boo, which has few intrinsic ab- 
sorption lines except for two strong hydrogen absorption 
lines : Pa7 (~ 10935 A) in Y band (order 51 and 52) 
and Paf3 (~ 12815A) in J band (order 44). Before the 
division, we removed these hydrogen lines by fitting with 
Voigt profile, and scaling the telluric absorption lines to 
correct the airmass difference (J2]). Figure [3] shows the 
flux calibrated spectra of images A and B. We used the 
continuum of the telluric standard star 10 Boo (5.67 mag 
in J-band from SIMBAE0) for flux calibration, assuming 
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Fig. 3. — Upper panel: the whole Y- and J-band spectra of images A (lower spectra) and B (upper spectra) of B1422+231. The observed 
spectra were telluric-corrected and smoothed with an 11 pixel (~ 3A) box car to clearly show the spectra with an effective spectral resolution 
of R~4,000. The bumps on the continuum are broad Mg II and Fe II emission lines from the QSO itself at z =3.628. The vertical bars 
show the expected positions of absorption lines associated with the z = 3.54 system: solid and dashed lines show the positions of the 
detected and non-detected lines, respectively. Lower panel: the telluric absorption spectrum in the same wavelength range. This is made 
by normalizing the spectrum of the standard star, 10 Boo, after removing two hydrogen absorption lines (Pa/3 and Pa7) with Voigt profile 
fitting. 
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Fig. 4. — Upper panel: the normalized spectrum of image A around Mg II doublet lines of the z = 3.54 system after the telluric corrections. 
Lower panel: the normalized spectrum of the standard star 10 Boo, showing the telluric absorption lines. The strong O2 absorption band 
at around 12690 A apparently affects the detection and profiles of components 1 and 2 of Mg II A2796. 



the effective temperature of 9480 K. 

To avoid the influence of many bumps due to Fe II and 
Mg II emission lines from the QSO itself on the con- 
tinuum fitting, the normalized absorption spectra were 
made by fitting the continuum in a velocity range of ±500 
km s _1 , which is narrower than the velocity widths of the 
emission lines (<~ 1000 km s -1 ), around the detected ab- 
sorption lines using splinc3 function. 

Heliocentric corrections of —21.89 km s _1 and —5.18 
km s _1 were applied to Y- and J-band spectra, respec- 
tively. By comparing the detected Mg II absorption lines 
with C II and Si II absorption lines for image A (RSB99), 
a residual offset of —3.4 km s^ 1 (about a half pixel) was 
found in our J-band spectrum. We shifted the absorption 
line spectrum by +3.4 km s _1 to match the wavelength 



of the C II /Si II absorption lines with those derived 
by RSB99 because their spectral resolution (i? =70,000, 
At> = 4.4 km s _1 ) is much better than ours (R =10,000, 
Av = 30 km s _1 ). This offset is likely to come from the 
slight change of the instrument setting that occurred be- 
tween the observation of targets and the acquisition of 
the Argon lamp spectrum. For y-band data, no absorp- 
tion lines that have obvious peaks as Mg II absorption 
lines of the z = 3.54 system were detected. Therefore, we 
did not shift the spectra after the heliocentric correction. 

4. Z = 3.54 SYSTEM 

4.1. Velocity Components 

In the final J-band spectrum, strong Mg II AA2796, 
2803 absorption lines with two velocity components at 



5 



TABLE 1 
Detected Mg II absorption lines. 



Image 


Comp 


W r (A2796)(A) 


M/ r (A2803)(A) 


Doublet ratio 


A 


1 


(0.020) 


0.018 


(1.1) 




2 


(0.074) 


0.052 


(1.4) 




3 


0.191 


0.124 


1.5 


B 


2 


(0.080) 


0.040 


(1.6) 




3 


0.111 


0.056 


2.0 



Note. — W r of the components that arc significantly affected by 
telluric absorption lines arc shown with parentheses. 



z = 3.54 are clearly detected for both images A and B 
(Figure [3). The corresponding Fe II A2383 (/ = 0.320) 
and A2600 (/ = 0.2394) absorption lines of this system 
are also detected in image A (Figure [3|), while other Fe II 
absorption lines, such as A2344 (/ = 0.1142) and A2587 
(/ = 0.0691), are not detected probably because of their 
small oscillator strengths. Hereafter, we will focus only 
on this z — 3.54 Mg II system: the search for other 
systems and their results will be discussed separately in 
S. Kondo et al. (in preparation). Unfortunately, the 
Mg II A2796 line is overlapping with the telluric O2 band 
as shown in Figure [4] and the systematic uncertainty may 
be left in the profile of Mg II A2796 even after the removal 
of the telluric absorption lines (§3.3). The impact of the 
telluric absorption lines will be discussed in the following 
subsection (§4.2) for each velocity component. 

This z = 3. 54 sy stem was first detected by 
iSongaila fe Cowid (|1996ft as an H I system with N(R I) 
= 2.3 x 10 16 cm~ 2 in their high-resolution optical spec- 
trum with Keck HIRES. Later RSB99 detected various 
metal absorption lines also with Keck HIRES in spatially- 
separated spectra of images A and C. Figure [5] shows the 
Mg II absorption lines of the z = 3.54 system for both 
images A and B from our J-band data, and C II, Si II, 
and O I absorption lines of images A and C from RSB99's 
optical data. 

For the Mg II absorption lines, we clearly detected 
two velocity components for both images A and B at 
~ —80 km s _1 and ~ +90 km s _1 , which are identi- 
fied as "component 2" and "component 3" in RSB99. 
The Mg II absorption line that corresponds to "compo- 
nent 1" is also detected at ~ —130 km s^ 1 , but only for 
image A (Figure [5|) . The rest-frame equivalent widths 
(W r ) of detected Mg II AA2796, 2803 absorption lines 
are summarized in Table [TJ Note that the equivalent 
widths of lines that are affected by the strong telluric 
absorption lines arc shown with parentheses. The total 
rest-frame equivalent width of A2796 for this system is 
calculated as 0.28A. Therefore, this system is classified 
as a "weak Mg II system" which is defined as a system 
with VF r (2796) < 0.3 A (jChurchill et al.|[l999l) . The cor- 
responding Fe II A2383, A2600 absorption lines are also 
detected for component 3 but only for image A as shown 
in Figure [6] (left panel). The Mg I absorption line is not 
detected for any component for both images within the 
uncertainty. 

We fit a Voigt profile to the Mg II and Fe II absorption 
lines and measured the column density, th e relative veloc- 
ity, and the Doppler width using VPFITFl fCarswell et al.1 

5 VPGUESS is a graphical interface to VPFIT written by Jochen 
Liske; see http://www.eso.org/~jliske/vpguess 
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Fig. 5. — Spectra of the z = 3.54 absorption system of images A 
(black lines), B (gray lines for Mg II and Fe II), and C (gray lines 
for Si II, C II and O I) of B1422+231. The spectral resolutions are 
about 10,000 and 70,000 for our SUBARU data (top three panels) 
and RSB99's Keck data (bottom three panels), respectively. The 
vertical dashed lines show the peak position of the three velocity 
components of image A. The velocity is relative to z = 3.53850. 
The atomic species and their rest-frame wavelengths are shown at 
the left bottom corner of each plot. 

[19871) and VPGUESS^ software assuming each veloc- 
ity component is composed of just one velocity sub- 
component, which is a good approximation for estimat- 
ing column densities from our data with the velocity res- 
olution of Av = 30km s _1 at most. The results are 
summarized in Table [5] The shown uncertainty is only 
that for the VPFIT fitting. For non-detected absorption 
lines, we calculated the 3c upper limit of the column den- 
sity from the S/N of the continuum around the absorp- 
tion lines. We discuss the characteristics of the detected 
Mg II and Fe II absorption lines in the following. 

4.2. Mg II Absorption Lines 



6 VPFIT is a Voigt profile fitting package provided by Robert 
F. Carswell; see http://www.ast.cam.ac.uk/~rfc/vpfit.ht mll 
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Table [T] summarizes the Mg II doublet ratio, which is 
defined as W r (2796)/ W r (2803), for each component of 
each image. Where both AA2796, 2803 absorption lines 
are not saturated, the doublet ratio should be equal to 
2 by simply reflecting the oscillator strengths. For the 
strongest component (component 3), the doublet ratio 
is found to be 1.5 for image A, while equal to 2 for im- 
age B. Therefore, the absorption line of image A may be 
slightly saturated while that of image B is not. Although 
the other weaker components are also unlikely to be sat- 
urated in view of the smaller W / r (2803) than that for 
the component 3 of image B, the doublet ratios for those 
components show values less than 2. This is probably be- 
cause W r (2796) of components 1 and 2 is underestimated 
due to the incomplete removal of the heavy telluric ab- 
sorption lines that are overlapped on those components 
(Figured]). Here, we discuss details of the Mg II absorp- 
tion lines for each component 1,2, and 3 . 
Component 1. For this component, RSB99 shows that 
there is no difference between the absorption lines of im- 
ages A and C, thus the same absorption profile was ex- 
pected for image B, which is located in between images 
A and C on the sky. However, we could not detect com- 
ponent 1 in image B despite its existence in the spectrum 
of image A (Figure [5]). In particular, we had expected to 
detect the Mg II A2803 absorption line because it is not 
affected by the strong O2 telluric absorption lines unlike 
the Mg II A2796 absorption line. The Mg II A2803 ab- 
sorption line for component 1 is located right in between 
two weak telluric absorptoin lines (Figure |4]) and should 
not be affected by the removal of the telluric absorp- 
tion lines. Therefore, we conclude that the gas cloud of 
component 1 covers only images A and C on the sky, sug- 
gesting a small-scale sub-structure on this high-redshift 
cloud. 

Component 2. Although we tried to fit both Mg II 
AA2796, 2803 lines of component 2 simultaneously, the 
uncertainties of the three fitting parameters (column den- 
sity, redshift and Doppler width) were found to be very 
large and we could not get any reliable estimate. This 
is most likely due to the strong telluric O2 band on the 
Mg II A2796 line (Figure 2}: systematic uncertainties 
may be left in the profile of Mg II A2796 even after the 
removal of telluric absorption lines (§3-3), which is con- 
sistent with the strange doublet ratios for component 
2 described above. Therefore, we used only the Mg II 
A2803 line for the Voigt profile fitting for both images A 
and B for this component. The results are summarized 
in Table [5J Due to the lack of the information on Mg II 
A2796, the uncertainties are fairly large especially for the 
Doppler width. 

For image A, fortunately, we have the information on 
the Doppler width with a high precision from the high- 
resolution optical spectrum (RSB99). We tried to esti- 
mate the Doppler width of Mg II A2803 line from those 
of C II (9.8 ± 0.5 km s" 1 ) and Si II (6.7 ± 0.5 km s" 1 ) 
lines. Generally, a Doppler width b consists of thermal 
ingredient 6th and turbulence ingredient &tu' 

b 2 = b 2 th + b 2 tu (2) 

and the thermal ingredient is expressed as : 



where k is the Boltzman constant, T is the equilibrium 
temperature of a gas cloud, and m is the particle mass. 
From the Doppler widths of two species of different mass 
(carbon and silicon) , the temperature and the turbulence 
Doppler width were estimated as T = 65, 000±18, 000 K, 
and 6tu = 2.6 ± 2.7 km s _1 , respectively, resulting in the 
Doppler width of Mg II of 7.2 ± 1.3 km s" 1 . This value is 
consistent with the Doppler width estimated solely from 
the Mg II A2803 line (5±29 km s _1 ), but with much 
improved uncertainty. 

Component 3. Because the influence of the telluric 
absorption lines is little for component 3 (Figure SJ, 
which is consistent with the doublet ratio of compo- 
nent 3 for image B as discussed above, we fitted a 
Voigt profile to this component for both Mg II AA2796, 
2803 lines simultaneously (Table [2]). Because the opti- 
cal absorption profiles are asymmetric for both images A 
and C (see Figure [5]) , RSB99 fitted three velocity sub- 
components to the profile of image A. Because they did 
not show the fitting results in their paper, we newly fit- 
ted the three sub-components to their dataQ (Figure [7]) 
: the fitting results for the three sub-components are 
logiV(C II)[cm- 2 ] = 13.06, 13.49, 12.96, and b=3.5, 9.4, 
14.9 km s _1 , respectively. The observed Mg II absorp- 
tion profile of component 3 is also found to be slightly 
asymmetric. Although, at first, we tried to fit multi- 
components to component 3, Mg II absorption lines with 
the three components identified by RSB99, the resultant 
fitting uncertainties were very large because of the lower 
spectral resolution (Av ~ 30 km s _1 ) compared with the 
widths of three components (b < 10 km s _1 ). Therefore, 
we treat component 3 as a single component for both 
images in the fitting with VPFIT. Although the resul- 
tant Doppler width does not have any physical meaning, 
the estimated column density is expected to be accurate, 
at least for image B, because the Mg II absorption lines 
are not saturated for image B as expected from the dou- 
blet ratio. For image A, though the Mg II A2803 line is 
expected to be slightly saturated in view of the column 
density (10 13 cm~ 2 ) and the expected smallest Doppler 
width (9 km s _1 ), the resultant column density of Mg II 
is expected to be reasonably accurate. 

4.3. Fe II Absorption Lines 

Because iron is essential to assess the chemical abun- 
dance of the object, the detection of Fe II absorption 
lines for the z = 3.54 system is important to understand 
the nature of this high- z gas cloud. Note that the de- 
tection of Fe II A1608 in Keck HIRES spectra is not 
reported so far because the past optical observations o f 
B1422+231 (IRauch et al.l[2^01^[Song"aila fe Cowid[l996l) 
do not cover the wavelength of the shifted Fe II A1608 
absorption line of z = 3.54 system (~ 7300 A). In our 
Subaru IRCS spectra, two Fe II lines, A2383, A2600, with 
the largest oscillator strengths were detected, but only 
for component 3 of image A, which has a higher column 
density for all the ct-elements. 

In order to examine the other velocity components of 
the Fe II absorption lines for each image A and B, we 
combined spectra for A2344, A2383 and A2600 lines to 
improve the S/N for Fe II detection. We did not use 



(3) 



7 The machine-readable data were kindly provided by Dr. 
Rauch. 



7 



1.2 

1 

0.8 
0.6 
0.4 
1 

0.5 
1 

0.5 
1 

0.5 



Fell comb 
1 



Fell 2344 

■ 1 



Fell 2383 
— I 



w 



Fell 2600 



l3£ 




f =0.1 142 

1 ■ 



f=0.320 
h 



f=0.2394 



1.2 
1 

0.8 
0.6 
0.4 
1 

0.5 
1 

0.5 
1 

0.5 



I 

riJ L ^^~ ,L i_r 

Tell comb 


i 


i 


1 

Fell 2344 

1 1 




i 


1 

f=0.1142 

1 1 


1 ^HpnJ ST 

Fell 2383 
'- — 1 




1 


f=0.320 
1 1 


Fell 2600 

i 








i 


f=0.2394 

i 



0.5 





1 

0.5 


i i i 


Mgll2803j 1 
i 


2 Ti 
i 


A 

i 


Mgll 2803 1 

i 


2 3 

i 


B 

i 


-200 




Velocity [km/s] 


200 




-200 




Velocity [km/s] 


200 



Fig. 6. — Spectra of Fe II lines of the z = 3.54 system for images A and B. The spectrum of Mg II A2803 line is also shown at the 
bottom as a reference. The top spectrum "Fe II comb" shows the combined spectra of Fe II A2344, A2383 and A2600. The spe cies and the 
rest-frame wavelength are shown at the bottom left of each spectrum. The oscillator strength of each line from Morton (1991) is shown at 
the bottom right of each plot. The spectral resolution of Fe II A 2600 and Mg II A2803 spectra is R ~ 10, 000, while that of Fe II A2383, 
2344 and "Fe II comb" is R ~ 5000. The vertical dashed lines show the center positions of three velocity components of each images. The 
gray lines show the Fe II absorption profiles predicted from the Mg II profiles (see the detail in Section 4.3). 



TABLE 2 

The Results of Voigt Profile Fitting with VPFIT 



Species 


Component 


logAf A (cm 2 ) 


logAT s (cm 2 ) 


va (km s 1 ) 


ug(km s 1 ) 


fcjlfkm s ) 


6s (km s 1 ) 


Mg II 


1 


11.9±1.9 




-127±16 




2±55 






2 


12.4±0.1 


12.20±0.07 


-78±5 


-85±3 


7.2±1.3 a 


10 ± 9 






(12.5± 0.4) 




(-77±6) 




(5±29)~ 






3 


13.03±0.07 


12.66±0.04 


90.6±0.6 


102.9±0.6 


9±1 


7±1 


Fe II 


2 


< 12.1 


< 12.1 












3 


12.72±0.05 


< 12.1 


83±3 




23 ±6 




Mgl 


2 


< 11.1 


< 11.1 











Note. — The column densities, relative velocities from z —3.53850, and Dopplcr widths measured with VPFIT arc shown 
for both line of sight A and B. 

a Estimated from the Doppler width of C II and Si II. See the detail in Section 4.1 



A2587 for this combining because of the intrinsic weak- 
ness and the resultant low S/N. First we transformed the 
profiles of A2344, A2383 lines to that of A2600 using the 
following equation: 



s'(A) = l-fl-s(Ax )) 

\ A2344 or 2383 / 



2600 



2344 or 2383 



(4) 



where s(A) is the observed spectrum, s'(A) is the trans- 
formed spectrum, and / is the oscillator strength of the 
lines. Before the combining, the spectrum for Fe II 
A2600 (R =10,000, Aw = 30 km s _1 ) is smoothed to 
match the spectral resolution to that of A2344 and A2383 
(R =5,000, Av = 60 km s _1 ). Finally each spectrum is 



combined with weighting by square of the S/N of the 
continuum around the absorption lines. The resultant 
spectra are shown at the top of Figure[6] We also plotted 
the expected absorption lines with gray lines in Figure 
[6] assuming log7V(Mg II)/7V(Fe II) = 0.31, which is the 
value for component 3 of image A in Table [2] when we 
assume that the absorption line is composed of a single 
velocity component. As a result of Voigt profile fitting 
for these detected Fe II absorption lines, they are found 
to have a very large Doppler width (23 ±6 km s _1 , Table 

Li. 

For image A, component 3 is clearly detected again 
on the combined spectrum, confirming the detection. A 
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Fig. 7. — Upper panel: the Mg II A2803 absorption profile of 
component 3 of image A (R ~10,000). The gray line shows the 
fitted Voigt profile. Although this velocity component is known 
to be composed of three sub- velocity components (RSB99) , we fit- 
ted just one component because the Mg II absorption profile is 
convolved with the instrumental function of ~30 kms -1 which is 
much larger than the intrinsic width (~10 kms -1 ) of the velocity 
sub-components. Lower panehthe C II A1335 absorption profile 
of component 3 of image A (R ~ 70, 000). The black line shows 
the observed normalized spectrum from RSB99. The three black 
dashed lines show the fitted Voigt profiles while the gray line shows 
the combined profile. 

weak absorption line is newly noticed between the wave- 
lengths of component 1 and component 2. However, be- 
cause the center velocity has a large offset from that of 
the Mg II line, it is hard to confirm the detection at more 
than a 3c level with the present data. Similarly, for im- 
age B, two dips near the center velocities of components 
1 and 2 are found in the combined spectrum, but the 
detection is tentative because the expected absorption is 
much weaker in view of the corresponding weak Mg II 
absorption lines. 

For component 3 of image B, compared to the expected 
absorption profile shown with gray line in the right panel 
of Figure El the observed combined spectrum is almost 
flat and the absorption appears to be significantly less 
than the expected amount. The 3c upper limit is calcu- 
lated as logiV = 12.1 for this non-detected component. 
From this, the column density ratio Mg II to Fe II for the 
component 3 of image B is estimated to be > 0.5, which 
appears to be much larger than that of image A, 0.31. 
This may suggest a large variance of the column density 
ratio, log7V(Mg II)/JV(Fe II), even in a small scale of 
~ 8pc, which is the projected separation between sight- 
lines A and B at z — 3.54. 

5. DISCUSSION 

5.1. Small-scale Structure of Absorbing Gas Cloud at 

z = 3.54 

By comparing the results for images A and B, consid- 
erable differences of the column densities of Mg II ab- 
sorption lines (log(N A /N B ) ~ 0.20±0.12 dex and 0.37 
±0.08 dex for component 2 and component 3, respec- 
tively) and considerable velocity shears (\va — v B \ ~ 
7.0±5.8 km s -1 and 12.3±0.8 km s" 1 at z = 3.53850 
for component 2 and component 3, respectively) are 
found as shown in Figure [5] Considering the align- 
ment of three images A, B, and C, the differences were 
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Fig. 8. — Velocity profiles of the z = 3.54 system for images 
A(Mg II, R =10,000), B(Mg II, R =10,000), and C(C II from 
RSB99, R =70,000) from bottom to top. The velocity is relative 
to z = 3.53850. The vertical dashed lines show the peak velocity of 
components 2 and 3 of image A. The horizontal dashed lines show 
the normalized continuum level. 

as expected from RSB99's interpretation of their opti- 
cal spectra of images A and C in that the relations of 
the column densities and the relative velocities among 
three images are log Na > log Nb > log Nc (assuming 
logiV(Mg II)/AT(C II) is equal for both images B and C 
) and \va\ < \vb\ < \vo\ for both components 2 and 3. 
The projected distances at z — 3.54 between images are 
<1ab = BAh^Q pc and cIac — 22. 2h^ pc. This very high 
spatial resolution (10 pc at z = 3.54 corresponds to ~ 
1 mas angular resolution for direct imaging) shows the 
power of the gravitational lensing (RSB99). 

B1422+231 also has a QSO absorption system at 
slightly higher redshift (z=3.624) and the transverse dis- 
tance reaches 1 pc, which is the smallest separation 
for QSO absorption systems ever studied with gravita- 
tionall y lensed QSOs (jBechtold fe Yedll995tlRauch et al.l 
l2001al) . though the system shows few variations of ab- 
sorption lines amon g multiple images (see Figure 5 of 
iRauch et al.l l2001al ) and this system is likely to be as- 
sociated with the QSO itself. Therefore, the ~ 10 pc 
structure for the z = 3.54 system is the smallest spatial 
st ructure ever detecte d for QSO absorption systems. 

Ellis on et alj (|2004l ) observed three lensed images of 
gravitationally lensed QSO APM08279+5255(z=3.911) 
with Hubble Space Telescope (HST). They detected many 
metal absorption lines at 1.1 < z < 3.8, which corre- 
spond to the transverse distance, from 30 Hjq pc to 2.7 
h^Q kpc. They found large variations of EWs for lower 
ionization systems, which are traced with Mg II dou- 
blet lines, even on the spatial scale of a few hundred 
pc while the higher ionization systems, which are traced 
with C IV doublet lines, show less va riations of EWs 
(jRauch et alj|2001at lEUTson et al.ll200l . Therefore, low 
ionization systems should reflect small scale gas clouds, 
which arc likely to be related to star formation activities 
in galaxies. Because this z — 3.54 system's spatial scale 
(~ 10 pc) corresponds to that of the typical Galactic 
gas clouds, such as giant mole cular clouds, H II regions, 
and SNRs (jRauch et al.lll999f) . the z = 3.54 system to- 
ward B1422+231 is a very precious target that enables 
the study of such Galactic-scale objects in detail at the 
galaxy-forming epoch. 

5.2. Expanding SNR Shell 
5.2.1. Shell Model 
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Fig. 9. — Images of the expanding shell to illustrate the formulation in the text (Section 5.2.1) for the case where information from three 
sightlincs, A, B, and C, is available. Left panel: the 3D image of sightlines that pass through the shell. Middle panel: the expanding shell 
seen from the direction of sightlines, showing the definitions of dj\, dg, dc and cj>. Right panel: the cross-section of the shell that contains 
the sightline of image A and the center of the shell, showing the definition of 6 a. Similarly, 8g and Qq are defined for the sightlines of 
images B and C, respectively. 
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Fig. 10. — Images of the expanding shell to illustrate the formulation in the text (Section 5.2.1) for the case where information from two 
sightlines, A and C, is available. Left panel: the 3D image of sightlines that pass through a circle that is off the center of the shell. Middle 
panel: the expanding shell seen from the direction of sightlines, showing the definition of 0. The radius of the cross-section is i?sin0. 
Right panel: the cross-section of the shell that contains the sightlines of images A and C. Though we can constrain i?sin© and usin©, we 
cannot infer the radius R and expansion velocity v due to the lack of information of another independent sightline, such as the sightline B. 



The column densities and the velocities of Mg II ab- 
sorption lines of components 2 and 3 vary among three 
sightlines as logiV^ > logiVe > log Ac and \va\ < 
\ v b\ < \ v c\- What types of object can generate these 
systematic variances on such a small scale? Based on 
the nearly symmetric velocity profiles of components 2 
and 3 with respect to the systemic velocity that corre- 
sponds to z — 3.53850 that are seen in optical spectra of 
both images A and B (see the profiles for C II and Si II 
in Figure [5]), RSB99 suggested that the components 2 
and 3 of the z — 3.54 system can be an expanding shell, 
such as an SNR. When the sightline passes through the 
center of a gas cloud of the expanding shell, the observed 
velocity becomes highest because the gas expands along 
the direction of the sightline while the column density 
becomes lowest because the path of the sightline in the 
shell is shortest at the center. On the other hand, when 
the sightline passes through the outer edge, the observed 
velocity becomes lowest, while the column density be- 
comes highest. Applying this model to the z = 3.54 
system, RSB99 found that the sightlines of images A 



and C pass near the outer edge and near the center of 
an expanding gaseous shell, respectively. The absorption 
lines of image B, which is newly observed in this study, 
shows intermediate values for both column density and 
relative velocity; this systematic kinematics seen in this 
gas cloud support the idea of RSB99 that this z = 3.54 
system is truly an expanding shell. Although a contract- 
ing (or collapsing) shell is an alternative choice, such an 
astronomical object is unlikely because the central ob- 
ject of the shell has to pull all the gas symmetrically in a 
subtle manner: normally such gas should fall through a 
disk or infalling envelope (not a shell). In fact, no such 
object is known in the Galaxy and near-by galaxies (see 
the listed exam ples o f astronomical objects in Section 5.2 
in lRauch et"alll2002f) . 

Following RSB99 and lRauch et~aH (|2002l) , we will con- 
sider a model of a three-dimensional (3D) expanding shell 
(Figured]) that has a radius of R and an expansion veloc- 
ity of i>exp- Those parameters can be constrained by two 
kinds of observables: the physical distance between two 
sightlines at the redshift (e.g., cLab) and the velocity dif- 
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ference of the two absorption components seen in a single 
sightline (e.g., Ava)- In case only two sightlines are avail- 
able, only the parameters for the expansion "ring" that 
is on the plane of the two sightlines can be constrained: 
one is the radius of the ring (i?sinO) and the other is 
the expansion velocity of the ring (i> exp smO), where O 
is the declination of the ring on the sphere (Figure [TU]). 
Although the relation between R sin and v cxp sin can 
be obtained, the 3D radius (i?) and the expansion veloc- 
ity (uexp) can never be determined because of the com- 
plete lack of the information on the absolute location of 
the ring (0) on the sphere: only the possible range of 
R sin and the lower- limit of u exp sin can be obtained 
(RSB99). 

On the other hand, in case three independent sight- 
lines (not on a single plane) are available, we can put a 
strong constraint on the expanding sphere. Since the ab- 
solute location of the planes that contain the sightlines 
cannot be determined in a unique way, R and v exp still 
cannot be determined. However, if the R value is fixed, 
the other parameter v cxp is determined because two 0s 
for the sets of two sightlines (e.g., A/B and B/C) can be 
determined from the two equations for two rings. There- 
fore, the relation between R and u CX p, which is useful 
to elucidate the astronomical nature of the shell, can be 
obtained as a fi nal product (see th e detaile d description 
in Ap pendix of IRauch et aLll2002D . In fact. IRauch et all 
(2002) suggested that the z = 0.5656 absorption sys- 
tem in the three sightlines of gravitationally lensed QSO 
Q2237+0305 is also an expanding shell because the ab- 
sorption lines showed variances similar to those of the 
z = 3.54 system. They analysed the relation between 
the radius and expansion velocity of the 3D shell and 
concluded that the expanding shell at z = 0.5656 can 
be interpreted as a supershell or a superbubble of 1-2 
kpc size that is frequently observed in the Galaxy and 
extra-galaxies. 

Based on the radial velocities in the two lines of sight 
(images A and C), RSB99 managed to constrain the 
parameters in a ring for the z = 3.54 system as 9 
h^Q < i?sinO < 34 h^g pc and u oxp sinO > 98 km s~ . 
Now, with the additional information of image B, we 
can obtain the relation of R and v cxp of the expanding 
shell at z = 3.54. Note that the model parameters can be 
completely determined if we have four independent sight- 
lines. Because B1422+231 has four independent gravita- 
tionally lensed images, future observations of the fourth 
image would be very valuable. 

We formulated the geometry of the expanding shell as 
shown in Figure [9] to obtain the radius- velocity relation 
of the z = 3.54 system. The difference of velocities of 
components 2 and 3 in the sightline of image A, (Av)a, 
and the distance from the center of the shell to the sight- 
line A, cLa, can be expressed as 

{Av) A = 2v cxp cos9 A (5) 
oIa = Rsin0A 7 (6) 

where 9a is the angle from the center of the shell to the 
sightline A (see Figure The equations for images B 
and C can be given similarly, resulting in six equations 
in total. Next, we can derive two equations about the 
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Fig. 11. — Relation of radius and expansion velocity of the ex- 
panding shell model. The solid line shows the relation for the 
z = 3.54 system in the case of (Av)g = 188 km s -1 , which is the 
most likely value from the observation (see Section 5.2.1 for de- 
tail). The dotted lines along the solid line are the case of (A»)g = 
185, 191 km s — 1 , showing la uncertainty of our calculation. The 
filled and open circles show the radius and expansi on velocity of 
the observed SNRs and H II regions in our Galaxy ||Kgo_fc_JIeiles] 
1991; Shields 1990; Kothcs & Kcrton 2002; Arnal & Corti 2007; 
Daiglc ct al. 2007; Cappa ct al. 2008; Cichowolski et al. 2008). 

geometry of triangle ABC : 

d\ = d B + A~B~ 2 + 2d B ABcos<j> (7) 

d 2 c =d 2 B + BC 2 + 2d B ~BC cos [ZABC - <f>], (8) 

where AB and BC are projected distances between im- 
ages at z — 3.54, cf> is the angle between BA and BO, 
and ZABC is the angle between BA and BC. Now, the 
parameters (Av) A , (Av) B , (Av) c , ~AB, BC, ZABC, are 
observed and the nine parameters R, v cxp , dA, d B , dc, 
9a, b , 9c, 4> are not determined. If these eight equa- 
tions are combined, the R(v cxp ) relation can be obtained 
as a solution^- 

For the z — 3.54 system, the projected distances are 
AB = 8.4/ifg 1 pc, BC = 14.4/ifo 1 pc, and AC = 22.2^ 
pc at z = 3.54. The angle between AB and BC is 153 . 
The velocity differences between components 2 and 3 are 
(Av)a = 168 ± 1 km s _1 , (Av) B = 188 ± 3 km s _1 , 
and (Av)c — 197 ± 1 km s _1 . By putting those val- 
ues into the equations, we obtained the function R(v exp ) 
that is shown with a solid line in Figure Ql] along with 
two dashed lines that show the la uncertainty range of 
our calculation : each line corresponds to the cases of 
(Av) B = 185 km s _1 , 191 km s _1 . Here, we are only 
concerned with the uncertainty of (Av) B because it is the 
largest among six observables : the available spectrum of 
image B is only our lower-resolution near-infrared data 
while the available spectra of the other two images A and 
C are higher-resolution optical data of RSB99. 

5.2.2. The SNR Origin of the Shell 

Figure [11] shows that the radius of the shell is larger 
than ~ 20 pc and the expansion velocity is larger than 
~ 120 km s" 1 . Although there is no a priori constraints, 

8 T he analytic solution is given in Equation (A5) in IRauch et al.l 
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the shortest radius (R < 30 pc) is unlikely because it re- 
quires: (1) an exact alignment of the shell with the three 
sightlines that have similar separations as the shell diam- 



ity are modeled as 



eter and (2) very high velocity (w exp > 200 km s^ 1 ) that R ^ ~ 21 - 6 P C ( i()5ierg ) \ 1cm" 3 



appears to be too fast for an H I shell. The largest radius 
(R 3> 100 pc) is also unlikely since it requires similarly 
very high expansion velocity and all three sightlines must 
pass near the edge of expanding shell, exactly. Therefore, 
the bottom of the curve (R ~ 50 — 100 pc, v cxp ~ 130 
km s _1 ) would be the natural choice of the radius and 
the expan sion velocity of th e shell (see the original dis- 
cussion in lRauch et al.ll2~002l for another QSO absorption 
system) . 

To compare to the known astronomical shell-like ob- 
jects, we also plotted the radii and expansion veloc 
ities o bserved for the Galactic H II regions ( Shield; 
19901 iKothes fc Kertonl [2002] lArnal fc Corti 



En 



no 



"45 / V exp (t) 



Dai gle et al J 120071 iCappa et al.fl200 8; Cichowols 



fields! 
2007T 



et al.1 



2008(1 and SNRs that have a struct ure of expanding shel l: 
the data points of SNR are from IKoo fc Heilesl (jl991l) . 
who observed 15 old SNRs with H I 21cm line in the disk 
of our Galaxy. As a result, the location of the z = 3.54 
system on the plot is found to be quite consistent with 
SNRs while all the H II regions show much smaller ex- 
pansion velocity less than 20 km s _1 . In fact, the distri- 
butions of SNRs match well with the bottom of the curve 
for the z = 3.54 system, which is quite consistent with 
the above consideration that the bottom of the curve is 
the most likely location of the shell. Therefore, we con- 
cluded that the shell is truly an SNR shell, confirming 
RSB99's initial suggestion in more solid way. 

5.2.3. Physical Parameters of the SNR Shell 

Assuming R = 50 — 100 pc and v cxp — 130 km s _1 
as the radius and the expansion velocity of the observed 
shell, we compare the observed SNR shell with the typ- 
ical self-similar expansion model to estimate the main 
physical parameters of the shell. Then, we compare the 
estimated physical parameters to those of typical SNRs 
to check the consistency of the SNR interpretation. In 
the following, we will utilize t he fo rmulation desc ribed in 
text-books by iTielensI ([20051 ) and iDraind (|2011[ ) to first 
estimate the age (t) of the SNR shell, and finally the total 
energy (E) of the shell, along with an independent esti- 
mate of mass (M) of the shell from the column densities 
of the absorption lines. 

If the self-similar concept is introduced with a simple 
power laws (Roct n ), the age of the shell can be roughly 
estimated from R and v cxp as t ~ 10 5 yr for z — 3.54 
system since t ~ R/v exp . The SNR shell of this age is in 
the radiative expansion phase (snowplow phase) rather 
than adiabatic expansion phase (Sedov phase). In this 
case, rj = 2/7 and the age t can be derived with 



f = 1.1 x 10 5 yr 



( R 



V50pcy V 130km s" 



(9) 



For the observed v exp (130 km s" 1 ) and the range of R 
(50 - 100 pc), the age of the z = 3.54 shell is estimated 
as: 

1.1 x 10 5 yr < t < 2.2 x 10 5 yr- (10) 
In this snowplow phase, the radius and expansion veloc- 



V 250 km s- 1 
. (U) 

where Eq is the total energy of a supernova, hq is the 
number density of the interstellar gas before the super- 
nova explosion. 

Next, we will discuss the mass of the shell (Af to t)- 
RSB99 estimated the number density (n) and the size 
along the line of sight (L) for the component 3 of image 
A using photoionization model as 



0.16 cm 3 < n < 1.6 cm 
0.015 pc < L < 1.6 pc. 



-3 



(12) 
(13) 



From these parameters, RSB99 estimated the mass of a 
gas cloud (M tot ) assuming two geometrical cases ; first 
is a homogeneous cylindrical slab with a thickness L and 
radius AC as a lower limit of the mass; second is a spheri- 
cal cloud with a radius AC as an upper limit of the mass: 



0.4M Q < Aftot < 2700M, 



0- 



(14) 



Here, we newly obtained an additional constraint on the 
radius R of the gas cloud as R ~ 50 — 100 pc. Assuming 
that the gas cloud has a shape of spherical shell with an 
average radius of R, thickness of L, and number density 
of n, we can calculate the total mass M to t of the shell 
with 



M tot = 47ri? Lnm H m, 



(15) 



where m H j is the mass of a hydrogen atom particle 
( m Hl — 1-66 x 10~ 24 g) and \x is the reduced mass 
(fi = 4/3). Here we ignored the effect of cos 9a on 
the thickness of the shell (see Figure [S]) in view of the 
large uncertainty of L. Using the value range of n and 
L (Equation [12] and [T3| and our result on the radius 
(50 < R < 100 pc), Mtot is estimated as: 



25M < M tot < 1050M, 



©• (16) 

This is consistent with the original mass estimate by 
RSB99 (Equation [T4]) , but narrows the mass range by 
two orders of magnitude. 

Because the uncertainty of this estimate is quite large, 
we try an alternative mass estimate in the following. We 
will estimate the mass in two ways here, using the column 
density of H I or Mg II, in order to check the consistency. 
In both estimates, we will use the column density of com- 
ponent 3 of image A because the ionization parameter, 
logf/, of this component was specifically estimated as 
log U = —4.4 (RSB99) which is low enough that we can 
estimate the mass simply from observed column density 
without any ionization correction. Although the column 
density of this component may not be the representative 
value of the shell, we would not expect a large uncer- 
tainty of more than an order of magnitude in view of the 
column density variation among the components seen in 
images A, B and C. 

First, the H I mass of the shell can be simply calculated 
as : 



M Hl = 47ri2 2 2V Hl m Hl 



(17) 



where R is the size of the shell and iV H j is the observed 
H I column density. Using log A^(H I)=16.05 ± 0.15 for 
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component 3 of image A (RSB99), the H I mass is calcu- 
lated as 

2.8M© < M Hl < 11M (18) 

Because this gas cloud is optically thin, almost all of the 
hydrogen is likely to be ionized. To estimate the total 
mass of the shell from the H I mass, we must eva luate 
the degree of ionization. In lDonahue fc Shulll ()199lD . the 
fraction of H I is calculated with 



n(H I) 

n H 



= (4.6 x IQ- 6 )U 



-1.026 



-4.7 < logt/ < -1. 



( 19 ) 

where n(H I) is the number density of only H I and njj is 
the total number density of hydrogen atom that includes 
H II. Since the ionization parameter of the component 
3 of image A is estimated as \ogU = —4.4 (RSB99), 
ri(H l)/n(H) is calculated as 0.15. Then, the total hy- 
drogen mass is estimated as: 

19M < M H < 75M . (20) 

Finally, the total mass of the shell can be calculated by 
multiplying the reduced mass, fi, as: 



25M < M tot < 99M, 



(21) 



This range is consistent with the typical scrambled gas 
mass of the observed SNR (10 - 1000 M^), with a r adius 
from about 10 to a few 100 pc (|Koo fe Heileslll991l ). 

Next, we attempt one more independent estimate of 
the total mass of the shell based on the column density 
of Mg II instead of H I. The total Mg II mass in the shell 
can be estimated using Equation (|17[) but after replacing 
H I to Mg II. For column density AT M jj, we use the 
value of component 3 of image A because this component 
is examined in detail by RSB99. As a result, the total 
Mg II mass is calculated as 



0.06M Q < M M 11 < 0.24M G 



(22) 



for the assumed range of radius. 

The total mass of the shell can be estimated first as- 
suming that all the magnesium is in the form of Mg II. 
This assumption is reasonable because the Mg I absorp- 
tion lines are not detected (the 3<r upper limit is calcu- 
lated as log N(Mg I) [cm' 2 ] < 11.1; see Section 4.1 or Ta- 
bled Although the possibility of the existence of a signif- 
icant amount of Mg III cannot be dismissed, we ignored 
the higher ionization states because the ionization pa- 
rameter log U of this component is estimated to be quite 
low as log U = — 4.4 (RSB99) based on th e photoioniza- 
tion modeling of iDonahue fc Shulll (|1991f l. We assumed 
the solar abundance, which is suggested by RSB99 for 
component 3 based on t he photoionization modeling by 
IDonahue fc Shulll (|1991| ). Wi th the solar abundan ce of 
magnesium (0.13% in mass; iGrevesse et"aT1l2010D . the 
total mass M to t is estimated from Equation (|22[) as 

47M Q < M tot < 188Af Q . (23) 

This mass range is pretty much consistent with the esti- 
mate from H I column density, Equation (I2ip . 

From the estimated radius and expansion velocity, we 
can finally constrain the energy of supernova explosion 
using Equation (|11|) . The remaining parameter in this 
equation, no, which is the number density of interstellar 
medium around the supernova before explosion, can be 



+ + + ++ ++ + + 

+ +VV D+ + 

+ + + + v . 
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Fig. 12.— JV(Mg II) vs. j V(Mg Il)/jV(Fe II) of we ak Mg II sys- 
tems. T he crosses ar e from Narayanan et al. (2008), the squares 
are from Rigbv ct al. (2002), and a filled circle shows the compo- 
nent 3 of image A of the 2 = 3.54 system. The filled squares show 
confirmed Fe-rich s ystems while a op en square shows confirmed 
non Fe-rich system jRigbv et al. 2002). The dotted line shows the 
solar abundance ratio (Asplund ct al. 2005). The 2 = 3.54 sys- 
tem (filled circle) is located in a horizontal branch near the solar 
value, where the confirmed Fe-rich systems (filled squares) are dis- 
tributed. 

estimated assuming that the shell consists of all of the gas 
that existed in the sphere with radius R before explosion, 
with the following equation: 



3M tl 



n 



ATrR 3 rriHfJ- 



(24) 



With the range of 50 < R < 100 pc and M tot from Equa- 
tion (|23l) . n is calculated as: 



1.8 x 10" 3 cm- 3 



< n < 3.6 x 10~ 3 cm~ 3 



(25) 

Then, the energy of supernova explosion, Eq, is calcu- 
lated with Equation ([TT]) : 



3.8 x 10 49 erg < E < 3.2 x 10 5U erg 



50, 



(26) 



The estimated energy is roughly consisten t with the en- 
ergy of supe rnova explosion, ~ 10 51 erg (jTielend 120051 : 
Drainc 2011)). The slight difference can be attributed to 
the l eft-over gas inside the shell (see e.g., IK00 fc Heilesl 
fl99lh that can effectively increase no, thus E through 
Equation (TTTjl . 

In this subsection, we have gone through the physical 
parameters of the z = 3.54 system as an expanding shell 
of an SNR. The expanding shell model and estimated 
physical parameters appear to be quite consistent with 
the properties of SNRs observed in the galaxy. With our 
calculation, this system is likely to be an SNR of about 
0.1 Myr with a radius of 50 — 100 pc, an expansion 
velocity of about 130 km s , and the total energy of 10 50 
erg. Therefore, we conclude that this z — 3.54 system is 
truly an SNR. 
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5.3. Type la Supernova ? 

Recall the very broad profile of Fe II absorption lines 
of component 3 in image A, which is described in Section 
4.2. What does this feature mean in the SNR interpre- 
tation of the z = 3.54 system? To answer this question, 
we first discuss the iron richness of the SNR shell, then 
examine the broad absorption feature in a more rigorous 
way to suggest that the iron is localized in the SNR shell. 
Finally, we conclude that the SNR shell is related to an 
SN la. 

5.3.1. Fe Richness 

The amount of iron in the gas cloud is crucial to in- 
vestigate the chemical enrichment by supernovae. Fig- 
ure [ft] shows the distribution of log iV(Mg Il)/AT(Fe II) 
versus logiV(Mg II) for weak Mg II sy stems. The 
crosses are from INaravanan et all (|2008f ) . who stud- 
ied 100 weak Mg II systems at 0.4 < z < 2.4 us- 
ing VLT data. The four squares are weak Mg II sys - 
tems studied with Keck data by iRigbv et al l {2002). 
Note that neither sample includes the systems whose 
Fe II are not detected. The filled and open squares 
are confirmed Fe-rich and non-Fe-rich systems based on 
their photoionization modeling using CLOUDY, respec- 
tively. The filled circle shows the z = 3.54 system 
(logiV(MgII)/Ar(FeII)=0.31±0.07). The dotted line 
show s the solar abundance ratio, log[7V(Mg)/7V(Fe)] Q = 
0.08 (jAsplund et al. II2005D for the case that all Mg and 
Fe atoms are in the Mg II and Fe II ionization states. 
The value of the z = 3.54 system is found to be closer 
to the solar value than those f or the other weak Mg I I 
systems with similar iV(Mg II). Narayan an et al.l (2008) 
suggested that any system near the solar value (dashed 
li ne in Figure HH ) arc truly Fe-rich systems. 

IRigbv et alT 72002) suggested that their three Fe-rich 
systems (iV(Fe II)~iV(Mg II)) have small sizes of ~ lOpc, 
and high metallicity of > 0.1 -Z© (see filled squares in 
Figure IT2")) . They specifically predicted that the z — 3.54 
system toward B1422+231 should show strong Fe II lines 
(iV(Fe II) ~ AT(Mg II)) in view of the small spatial struc- 
ture (^10 pc) and high metallicity inferred by RSB99. 
In fact, log Af(Mg Il)//V(Fe II) of the z = 3.54 system is 
found to be similarly Fe-rich as the three Fe-rich systems 
(see Figure [T^]) , suggesting solar to sub-solar metallicity. 

All those results support the iron-richness of the SNR 
shell at z = 3.54. The high iron abundance of the 
SNR naturally su ggests that it is an SN la. In fact, 
IRigbv et al.1 ([20 02) suggested that the three iron-rich sys- 
tems in their samples have [a/Fe]< and have been en- 
riched by SNe la because high iron column density that 
is similar to magnesium column density cannot be ex- 
plained by other enrichment processes such as SNe II. 
Therefore, it is highly likely that the z = 3.54 system is 
a gas cloud enriched by SNe la. More detailed arguments 
based on abundance estimate using CLOUDY photoion- 
ization modeling will be presented in our separate paper 
(S. Kondo et al. in preparation). 

If the gas cloud was truly enriched by SNIa explosion, 
the total amount of the iron in the gas should be consis- 
tent with the yield of the iron from SNIa explosion. As- 
suming the Fe II absorbing gas is in the form of shell with 
R = 50 — 100 pc as Equation (fTTj) . the Fe II mass is esti- 
mated as 0.07-0.29 Mq from the component 3 of image A. 
Our estimated Fe II mass is consistent with the observed 



mass of SNela. iScalzo et ail (|2010D and iSilverman et al.l 
(|201lD suggest that the estimated mass of radio active 
56 Ni (eventually decays to 56 Fe) ejected from an SN la 
ranges from 0.02 M to 1.7 M Q . Note that the slightly 
low estimated energy of the SN (Section 5.2.3) also fa- 
vors an SNIa interpretation rather than other types of 
SN with more energetic explosions. Because our esti- 
mate does not include the iron in other ionization states, 
such as Fe III, within the shell as well as the iron inside 
the shell, the total mass of iron is expected to be more 
than the estimated value. 

5.3.2. Fe Localization 

Another characteristic of the Fe II absorption line is 
the Doppler width (23 ± 6 km s^ 1 ) that appears to be 
unusually broad compared with that of a-element (e.g., 
^Mgll ~ 9 km s _1 ). This is quite strange for a QSO 
absorption system, because the mass of an iron atom is 
much heavier than that of a elements, thus the width of 
an Fe II absorption line should be smaller than that of 
a elements. We first compared the widths of Mg II and 
Fe II lines of the z = 3.54 system with the past surveys 
of weak Mg II systems (IRigbv et al.ll 2002HChur chill et al.l 
l2003bHNaravanan et al.ll2008fl . Figure [TBI shows correla- 
tion of the Doppler widths of Mg II and Fe II lines of 
weak Mg II systems. The points should be located below 
the dotted line (& Mg jj = b Fe jj) because iron is heavier 
than magnesium. In fact, most data points from the lit- 
eratures are distributed along or below the dotted line. 
However, the z = 3.54 system is found to be located sig- 
nificantly above the line even considering the uncertainty, 
suggesting the uniqueness of this system. 

We first checked the possibility of blending of other 
absorption lines on the Fe II absorption lines. For exam- 
ple, if Mg II absorption systems existed at z =3.22 and 
z =2.87, strong metal absorption lines Mg II AA2796, 
2803 would blend with the Fe II A2600 (A = 11804 
A) and A2383 (A = 10814A) lines, respectively. How- 
ever, such absorpt i on sys tems have not been repor ted in 
iSongaila fe Cowid (|1996t) and lRauch et al.l (|2001al) . who 
detected many Lyman series lines and C IV absorption 
lines of B1422+231. Therefore, we conclude that the 
broad features of Fe II lines are real. 

Here we try to evaluate the excess Fe II absorption by 
decoupling the broad feature into two sub-components. 
If we compare the Mg II and Fe II profiles, the excess ex- 
ists on the blue side of the Mg II absorption peak (Figure 
[6]). Then, we fitted two velocity components to both Fe II 
A2600 and A2383 lines with the fixed Doppler widths of 
9 km s _1 , which is the upper limit estimated from the 
Doppler width of Mg II (~ 9 km s _1 ). The results are 
shown in the bottom panel of Figure [14j The black line 
shows the component whose redshift is fixed to the value 
of the Mg II absorption line during the fitting, while the 
gray line shows the other excess component, for which 
the redshift was not fixed. The resultant difference of 
peak velocity between the two velocity components is 
26 ± 7 km s , and the column densities of blue and red 
components are log N [cm -2 ] = 12.3±0.2 and 12.5±0.1, 
respectively. This fitting appears to match well with the 
observed data for both Fe II A2600 and A2383 lines. 

Now, the question is what is the blue component ? 
Because there is no obviously corresponding a-element 
absorption lines, this iron gas cloud is inferred to be lo- 
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b Mgll [ kms ~ 



Fig. 13. — Correlation of Doppler widths of Mg II and Fe II 
abso rption lines of w e ak Mg II systems. The crosses are 
fromlRigbv et al.l J2002T) , IChurchill eTaLl J2003b) , Nar ayanan et al.1 
1120081) , while the filled circle shows the z = 3.54 system (compo- 
nent 3 of image A). The dotted line shows the location where the 
widths of Mg II and Fe II are identical. While the crosses are dis- 
tributed along or below the dotted line, the point of the z = 3.54 
system is located above the line with a large offset. 

calized in the shell. The localization of iron in the SNR 
also supports the SNIa origin of the SNR shell: the Fe II 
features imply the ejection of the mass from the SN la, 
while the other a-element absorption lines are likely to be 
the interstellar gas scrambled by the shock of the SNR. 
In our Galaxy, Hamilton et~aTl (|1997l ) observed SN1006, 
which is an SN la remnant, with absorption lines on a 
spectrum of Schweizer and Middleditch star, whose sight- 
line intersects near the center of SN1006. Despite the 
age difference (~1000 yr for SN1006, ~ 10 5 yr for the 
z = 3.54 system), it would be useful to compare those 
two objects in detail to infer the nature of the z — 3.54 
SNR shell because the geometrical configuration is very 
similar. For SN1006, very broad (~ 5,000 km s" 1 ) Fe II 
absorption line is detected at the systemic velocity while 
the ejected Si II shows a redshift of ~ 5,000 km s . This 
suggests that the mixing of ejected iron with interstellar 
gas occurs after the mixing of ejected a-elements. The 
two velocity sub-components of Fe II absorption lines of 
component 3 of image A can be interpreted as mean- 
ing that the red component with the same velocity as 
the a-element lines comes from Fe gas that mixed with 
the scrambled interstellar gas, and the blue component 
comes from Fe gas that has not been mixed with it yet 
(Figure[T5|). Therefore, we might witness the mixing pro- 
cess of ejected iron with interstellar gas at z = 3.54. 
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Fig. 14.— Voigt profile fitting to the Mg II A2803 (bottom), Fe II 
A2600 (middle) and A2383 (top) absorption lines of component 3 
of image A. The thin lines show the observed spectra. The solid 
thick lines show the fitted Voigt profiles with the Doppler width 
and redshift fixed to the values of Mg II absorption lines. The 
gray lines show the fitted Voigt profiles to the residual absorption. 
The dashed lines show the combined profiles of black and gray 
components. 
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5.3.3. SNe la at High Redshift 

In this section, we saw the evidence of the richness and 
localization of iron in the SNR shell at z = 3.54. The 
detected iron richness cannot be explained by processes 
other than SN la. Moreover, the suggested localization 
of iron in the shell supports the SN la interpretation. 
From these facts, we conclude that the z = 3.54 system 
is an SN la remnant. 

The extragalactic SN la has b een studied extensively 
for t he supernova cosmology (jGoobar fe Leibundgutl 
120111 ). However, even the most distant SN la event ever 
detected is at z = 1.55 (jRodnev et alJ[2012T ) and more 
distant objects that are important for the study of cos- 
mic chemical enrichment history are hard to detect and 
with 8-10 m class telescopes. For studying such high-z 



Fig. 15. — Schematic image of the SNR shell of the z = 3.54 
system. The Fe-rich gas clouds are localized in the expanding shell, 
while the a element gas clouds, such as Mg II, are distributed 
homogenously in the expanding shell. 

SNe la, absorption systems toward gravitational-lensed 
QSO may serve as good targets as is the case for this 
z = 3.54 system toward B1422+231. Even for the single 
sightline, the Fe-rich absorption systems as seen in Figure 
Q2] would become good targets for studying chemical en- 
richment history at high-redshift (z > 2) with more data 
available with sensit ive high-resolution spec troscopy with 
adaptive optics (see iKobavashi et"aT]|2005[ ). 

6. SUMMARY 
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We obtained near-infrared high-resolution (R =10,000) 
spectra of images A and B of gravitationally lensed QSO 
B1422+231 with the Subaru 8.2 m telescope and the 
IRCS echelle spectrograph. Although the observed PSFs 
of images A and B are partially overlapping in the slit, 
we managed to extract spectra for each image A and 
B to examine the differences of absorption lines. We 
detected Mg II AA2796, 2803 absorption lines of the 
z = 3.54 system, which had been found in images A 
and C in optical spectra by RSB99. Corresponding Fe II 
A2383, A2600 absorption lines are also detected but only 
for the component 3 of image A. The projected sep- 
aration between A and B images at z = 3.54 is just 
8.4/i-Tq 1 pc and we found considerable differences of col- 
umn density (A(logA) ~ 0.3 dex) and velocity shear 
(Av <~ 10 kms -1 ) between both images on such a small 
scale. These differences suggest the smallest structure of 
gas clouds ever detected for QSO absorption systems. 

Considering the physical origin of the differences of a- 
elements absorption lines among three images A, B, and 
C, we conclude that the z = 3.54 system is an expanding 
shell as originally suggested by RSB99. The information 
of three images enable us to analyse the 3D structure 
of the absorbing gas cloud and concluded that this z — 
3.54 system is an SNR shell whose radius and expansion 
velocity are 50—100 pc and <~ 130 km s _1 , respectively. 
We estimated several physical parameters of the shell, 
which are found to be consistent with SNRs: the age is 
~ 10 5 yr, the mass of the shell is ~ 100 Mq, and the 
energy of the SNe is ~ 10 50 erg. 

We also found that the Fe II absorption lines of 
z = 3.54 system have much larger column density 
and Doppler width than those of the weak Mg II sys- 
tems in the literature. The small column density ratio 
log iV(Mg II)/AT(Fe II) (=0.31±0.07) is indicative of the 
richness of the iron of the z = 3.54 system. We also 
estimated the mass of Fe II as 0.07—0.29 Mq which is 
roughly consistent with the yield of observed SNe la. 
Moreover, the large Doppler width of Fe II, which is inter- 
preted as the existence of one more velocity component 



that is blueshifted and extremely Fe-rich, suggests that 
the Fe-rich gas cloud is localized in the expanding shell 
of SNR. From these results, we conclude that the SNR 
shell at z — 3.54 is produced by an SN la explosion. 

Supernova explosions are thought to be one of the most 
important processes that drive the formation of galaxy 
through the dynamical energy input and the chemical 
enrichment. If this z — 3.54 system is truly an SNR, 
it is the farthest sample of a supernova ever observed. 
If higher S/N and/or higher spectral resolution data of 
B1422+231 were obtained in the future, the difference 
between iron and a-elements would be more clearly con- 
firmed. Then, physics of the mixing of interstellar gas 
with ejected gas from supernovae could be discussed in 
detail. If the spectrum of image D is also obtained, the 
radius and expansion velocity can be strictly determined 
with the expanding shell model. 

This z = 3.54 system illustrates the power of the 
gravitational lensing effect for the study of QSO absorp- 
tion systems with the extremely high spatial resolution 
and multiple sightlincs. In the case of present study, 
the separation between images A and B is only < 10 
pc at z=3.5 that corresponds to about 1 mas angular 
resolution. More observations of gravitationally lensed 
QSOs will reveal the kinematics (e.g., expansion) of gas 
clouds on interstellar scales at high redshift during cru- 
cial phases of galaxy formation. 
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